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Abstract 

The short-time behaviour of the critical dynamics for magnetic systems 
is investigated with Monte Carlo methods. Without losing the generality 
we consider the relaxation process for the two dimensional Ising and Potts 
model starting from an initial state with very high temperature and arbi- 
trary magnetization. We confirm the generalized scaling form and observe 
that the critical characteristic functions of the initial magnetization for the 
Ising and the Potts model are quite different. 
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1 Introduction 



The critical phenomena for magnetic systems in the equilibrium have long been 
studied. Due to the large correlation length near the critical point, the large scale 
behaviour of the systems does not depend on the microscopic details. An universal 
scaling form characterized by a number of critical exponents describes the scale 
invariance in critical systems. For magnetic systems in a non- equilibrium state, 
however, it is much less understood. Indeed, for dynamic systems in the long 
time regime where the equilibrium is almost reached, one finds also a dynamic 
scaling form. Due to critical slowing down, however, numerical study of critical 
dynamics is very difficult. Even for the simplest system as the two-dimensional 
Ising model, the dynamic exponent z has not rigorously been determined. Even 
though much effort has been put to develop fast algorithms to overcome critical 
slowing down, the dynamic scaling has not yet been deeply investigated. 

Recently a dramatic step forward in the study of critical dynamics has been 
made. For a relaxation process starting from an initial state with very high 
temperature and small magnetization, it is argued that universal scaling behaviour 
already emerges in the macroscopic early time |], @, §]• Important is that 
another dynamic exponent should be introduced to describe the dependence of 
the scaling behaviour on the initial magnetization. It has also been discussed 
that under certain conditions the effect of the initial magnetization remains even 
in the long-time regime of the critical relaxation [[5], ||. Numerical simulations 
support these predictions |7|, |9], [L(], [11], [12|, [L3|, [14], [T3| . In the equilibrium state, 
introducing the lattice size L as a scaling variable in the finite size scaling induces 
essential progress in the study of critical phenomana, especially in the numerical 
simulation. In dynamic systems, besides the lattice size L we now have the 
evolution time t as another scaling variable. Different from the case in the long 
time regime, it is demonstrated that in short-time dynamics the scaling variable 
t plays a very important role in understanding the scaling behaviour, especially 
in numerical simulations. The critical exponents can be obtained from the power 



law behaviour of the observables in the early time t [10, 12, [TI| 17] or from the 
time-dependent finite size scaling |TE|, [01. Numerical results indicate that the 
study of short-time dynamics not only enlarges the fundamental knowledge on 
critical phenomena but also provides new ways to determine all the static as well 
as the dynamic exponents. 

For the understanding of the universal behaviour in short-time dynamics, 
we should keep in mind that there exist two completely different time scales, 
the macroscopic time scale and the microscopic time scale. Universal behaviour 
emerges only after a time period which is long enough in the microscopic sense. 
This time period is called the microscopic time scale t mic . We expect that t m i C 
would be very small compared with the macroscopic time scale which is typi- 
cally t mac ~ L z or t mac ~ t~ uz . Here z is the dynamic exponent and r is the 
reduced tepemrature. However, how big t m i C is depends on the microscopic de- 
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tails. If t m i C would be so big as comparable with the macroscopic time scale t mac , 
what one claims for the universal behaviour in the short-time dynamics becomes 
meaningless. 



In Ref. [12 , 17 1, the short-time behaviour of the critical dynamics for the 



two dimensional Ising and Potts model has carefully and seriously been analysed. 
Simulations have been performed with both the heat-bath and the Metropolis 
algorithm. Fortunately, it is discovered that t mic ~ 5 — 30 Monte Carlo time 
steps for both algorithms. If one thinks a Monte Carlo time step is a typical 
microscopic time unit, this result is quite reasonable. Taking into account the 
possible effect of t m i C , critical exponents 6, z and (3/v are confidently obtained 



for both algorithms [|12| , [I7J . The results are well consistent and strongly support 
universality in the short-time dynamics. Interestingly, when the dynamic system 
exhibits already universal behaviour in the macroscopic short-time regime, the 
spatial correlation length is still not so big, e.g. for the Potts model at t = 100, 
£ ~ 5 [nj. Some more discussions on universality in short-time dynamics may 



also be found in the recent preprints |14|, |20|, |21], |22| . 

Does the critical exponent 9 (or xq) descirbe all the new physics in short-time 
dynamic systems? Very recently, it has been suggested that due to the fact that 
the time evolution of the magnetization is not a Markovian process, the global 
persistence exponent is also an independent critical exponent in the dynamic 
systems |23|, ^4[] . Numerical simulations show that this is indeed the case p3j , [20H . 
The power law decay of the global persistence probability is directly observed and 
the global persistence exponent for the two dimensional Ising and Potts model is 
determined to a satisfactory accuracy [ [20f . 

Since the scenario described above is so attractive and, on the other hand, 
all the results indicate that we are still far from a complete understanding of the 
critical short-time dynamics, we would like to push forward the study of more 
general dynamic processes. In a previous letter [|25|| , the author has numerically 
investigated the short-time behaviour of the critical relaxation for the two dimen- 
sional Potts model starting from an initial state with very high temperature and 
arbitrary magnetization. A generalized dynamic scaling is found to describe the 
crossover behaviour for the short-time dynamics from the repulsive fixed point 
mo = to the attractive fixed point mo = 1. The critical characteristic func- 
tion for the initial magnetization is numerically determined from the finite size 
scaling. 

The present paper is an extension of the previous letter. We aim to present 
a systematic description of the numerical data, and meanwhile pay attention to 
the microscopic time scale t mic . Compared with the data reported in the previous 
letter, the maximum evolution time t here is extended to t = 2000 for the lattice 
size L = 144 and t = 400 for the lattice size L = 72, twice as long as in the 



previous letter |25| . This allows us to see how stable the results are with respect 
to the time t. To reduce the fluctuations, we have also increased the statistics for 
the lattice size L = 144 for the bigger initial magnetization m . Furthermore, as 
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an alternative example, we have carried out also a similar simulation for the two 
dimensional Ising model. We interestingly observe that the critical characteristic 
functions for the initial magnetization for these two models are quite different. 
Finally some discussions are given corresponding to a recent preprint on this topic 

& * 

In the next section, we present scaling analyses for the short-time dynamics 
based on which the numerical data will be analysed. In Sec. 3, numerical results 
for the two dimensional Potts and Ising model will be reported. The last section 
contains the conclusions. 

2 Scaling analysis 

For the critical relaxation process of a dynamics of the Model A starting from an 
initial state with very high temperature and small magnetization, Janssen, Schaub 
and Schmittmann [0 proposed a dynamic scaling form which is valid up to the 
macroscopic short-time regime, 

M (k \t, t, m ) = b- kl3/ "M (k \b- z t, b 1/u T, b x °m ) (1) 

where M^ k ' is the /c-th moment of the magnetization, t is the dynamic evolution 
time, t is the reduced temperature, and the parameter b represents the spatial 
rescaling factor. Besides the well known static critical exponents (3, v and the 
dynamic exponent z, a new independent exponent xo which is the anomalous di- 
mension of the initial magnetization mo is introduced to describe the dependence 
of the scaling behaviour on the initial magnetizations. Based on the scaling rela- 
tion it is predicted that at the beginning of the time evolution the magnetization 
surprisingly undergoes a critical initial increase 

M(t)~m t e (2) 

where the exponent 6 is related to the exponent x by 9 = (x — f3/v)/z. 

The scaling form in Eq. (|1|) is valid only under the conditions that the ini- 
tial state is at very high temperature with small initial magnetization. For a 
critical relaxation from an initial state with very high temperature and arbitrary 
magnetization, is there universal scaling behaviour in the short-time regime? 

In Fig. [I], time-dependent magnetization profiles with different initial magne- 
tizations mo for the critical Potts model are displayed. The lattice size is L = 72. 
For smaller mo, the magnetization continously increases. The time scale for this 
incease is to ~ m z ^ x °. For bigger mo, the magnetization rises up rapidly and 
then decreases. After a macroscopic long time t 

mac ~ L z (out of the time range of 
Fig. |1|), we expect that all the profiles gradually join together and decrease expo- 
nentially. Traditionally, when we consider the dynamic scaling in the long-time 
regime the effect of the initial conditions is always neglected. However, as we can 
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see in Fig. [I], in the short-time regime the effect of the intial magnetization is 
very prominent, which is more or less due to the fact the exponent 9 is positive. 
Therefore it is important and interesting to study whether there exists universal 
behaviour in the short-time dynamics. 

A scaling form represents nothing else than a kind of scale invariance, or self- 
similarity in the systems. At the critical point the time correlation length of the 
dynamic systems is infinity, therefore we may expect a certain scale invariance 
independent of the initial conditions. The self-similarity of the magnetization 
profiles for different m may also be observed in In Fig. |I[ Important is that we 
should consider the scaling behaviour of the initial magnetization carefully. In 



the previous letter [25], the author proposed a generalized dynamic scaling form, 



e.g. for the fc-th moment of the magnetization, 

M {k) {t, r, L, mo) = b~ kp/v M {k) {b~% b 1/u r, b^L, <p(b, m )) (3) 

For the convenience of later discussions, finite systems have been considered here 
and L is the lattice size. The scaling behaviour of the initial magnetization mo is 
specified by the critical characteristic function ip(b, mo), which in the limit mo — > 
tends to the simple form b x °m used in Eq.(|l|), but is in general different. To avoid 
possible confusion, we have changed the notation x for the critical characteristic 
function in the previous letter to <p in this paper. This generalized scaling form 
describes the crossover behaviour between two fixed points, the repulsive one 
mo = and the attractive one mo = 1. It is also in a similar spirit as that in 
the correction to the scaling, where non-linear effects of an off-fixed point are 
considered. In general cases of the crossover phenomena, however, mo would also 
enter the scaling behaviour of the time t, r and the anomalous dimension of M^ k \ 
or the exponents z, f3 and v would depend on m . Then the situation will be 
very complicated and numerical simulations are difficult. 

The simple ansatz of the scaling form in Eq.(^) for the short-time dynamics 
that the exponents (3, v and z do not depend on the initial magnetization mo 
is based on the fact that if we believe there is a scaling form in the short-time 
regime of the dynamic relaxation process, this scaling form should smoothly cross 
over to that in the long-time regime. In other words, the initial magnetization 
should not enter the renormalization of the critical system in equilibrium and 
near equlibrium. However, the dimension of the initial magnetization would, 
in general, depend on the value of the initial magnetization. Assuming that 
we can either perturbatively or non-perturbatively perform the renormalization 
calculation of the dynamic system, the renormalization equation for the initial 
magnetization m may be written as 

d In m(A) , ,,,, , 

-a^t 1 =*(»>(*))■ w 

Here r] (m(\)) is the anomalous dimension of the initial magnetization and all 
the irrelevant operators are already set to the fixed points. The solution of the 
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equation (f|) gives the critical characteristic function ip(b, mo). Therefore <p(b, mo) 
can be written in a form 



^(b,m ) = F(bF-\m )) (5) 

with condition yj(l,m ) = m . If we define g = _F _1 (mo), under the renormal- 
ization group transformation g will be scaled as bg. Then we may conclude that 
the time scale where the effect of mo remains is g~ z . 

We should also point out that the dynamic scaling form in Eq.(]3|) may straight- 
forwardly be generalized to a critical relaxation process starting from an initial 
state with an arbitrary temperature and arbitrary magnetization. In general the 
initial temperature will also enter the dimension of the inital magnetization and 
vice versa. The critical characteristic functions for the initial magnetization and 
temperature are a bit more complicated. 

In order to verify the scaling form in Eq.(^), as a first approach, in this paper 
we always set r = 0. We perform simulations for a pair of lattice sizes L\ and L^- 
From the scaling collapse of the magnetization or its higher moments, in princi- 
ple, we can estimate the values of the function ip{b, mo) at b = L2/L1 for different 
mo- This works always well in case of smaller initial magnetization. When the 
initial magnetization m is bigger, the dynamic system leaves the repulsive fixed 
point m = and approaches to the attractive fixed point m = 1. The deter- 
mination of the function ip(b,m ) becomes more difficult since the dependence 
of the observables on the initial magnetization gets weaker. To reduce the extra 
fluctuations from M^ k \t, 1), we may introduce a magnetization difference 

M d (t, m ) = M« (t, 1) - M« (t, m ) (6) 

and a Binder-type cumulant 

M^(t,l)-M^(t,m ) 

^' m0) - [M d (t,m )P (7) 

to measure the function ip(b, m ). 



3 Generalized dynamic scaling 

In this section, we present numerical results for the short-time dynamics of the 
two-dimensional 3-state Potts and Ising model. In the first subsection, we system- 
atically analyse the data for the Potts model and demonstrate how to determine 
the critical characteristic function 93(6, mo) from the finite size scaling. In the sec- 
ond subsection, we report the results for the Ising model and meanwhile discuss 
some relevant important problems. 

In this section, all the simulations and discussions are restricted to the case 
r = 0, i.e. the dynamic system is exactly at the critical point. Therefore the 
exponent v will not enter the calculation. 
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3.1 The Potts model 



The Hamiltonian for the g-state Potts model is given by 

H = KJ2^ * = 1,...,<Z (8) 

<ij> 

where < ij > represents nearest neighbours. In our notation the inverse temper- 
ature is already absorbed into the coupling K. It is known that the critical points 
locate at K c = log(l + y/q). In this paper the 3-state case will be investigated. 
The exact value of the static exponents (3/v = 2/15 and the dynamic exponent 
z = 2.196(8) obtained from the power law decay of the auto-correlations ITO, [T2 



will be taken as input. In the numerical simulations, we measure the kth moment 
defined as 

M (k) (t,m )=<[^Y,(^m-l)] k >, (9) 

where the average is taken over independent initial configurations and the random 
forces. The initial magnetization m has sharply been prepared. 

In order to determine the critical characteristic function (p(b,mo), in the pre- 
vious letter we have performed the simulation for different lattice sizes, L = 
36, 72 and 144 with samples of different initial configurations 80 000, 40 000 
and 8 000 respectively. For the initial magnetization mo, we take the value 
m = 0.14, 0.22, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80 and 1.00 for the lattice size 
L = 144. For m > 0.50, the data are more fluctuating and therefore in this 
paper we increase the statistics for the lattice size L = 144 to 16 000 and for the 
lattice size L = 72 to 80 000. We updated the dynamic system to the maximum 
evolution time t = 2000 for the lattice size L = 144, t = 500 for the lattice size 
L = 72 and t = 100 for the lattice size L = 36. From the scaling collapse of an 
observable for two different lattice sizes L\ and L 2 with suitable initial magne- 
tizations, we can estimate the values of the function (p(b,mo) at b = L2/L1 for 
different mo- In principle, any observables which essentially depend on the initial 
magnetization may be used for this purpose. In this paper, the statistical errors 
are estimated by dividing the samples with different initial configurations into 
four groups. 



In Fig. |2a], the scaling plot for the magnetization is displayed for the lattice 
Li = 72 and the lattice L2 = 144. The initial magnetization for the lattice 
L2 = 144 is mo2 = 0.14. If one can find an initial magnetization moi for lattice 
L\ = 72 such that the time evolution curves of the magnetizations for both 
lattice sizes collapse, the scaling is valid and from the scaling form in Eq.(H) 



one gets (p(2, m^) = moi. In Fig. stars represent the time evolution of the 
magnetization of the lattice L2 = 144. Crosses are the same data but rescaled in 
time and multiplied by an overall factor 2^/". For the lattice L\ = 72, we have 
performed the simulation with two different initial magnetizations m = 0.175 
and m = 0.185, for which the magnetizations have been plotted by the lower 
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and upper solid lines. By linear extrapolation, we obtain the time evolution of 
the magnetization with initial values between these two values. The solid line 
laying on the crosses in Fig. 2a is the time evolution of the magnetization for 
L = 72 with an initial magnetization moi = 0.1802(6) which has the best fit to 
the crosses. Therefore, </?(2, 0.14) = 0.1802(6). Compared with the scaling plot 
for the magnetization difference Md(t,mo) presented in Fig. 1 in the previous 
letter, the evolution time has been extended to t = 400 in the time scale of the 
lattice size L = 72. The quality of the fit in Fig. [2a] is as good as that in Fig. 1 
in the previous letter. The resulting value for y?(2,0.14) = 0.1802(6) is also very 
well consistent with (^(2,0.14) = 0.1800(4) in the previous letter fl2*5| . 

However, when the initial magnetization mo is getting bigger, the determina- 
tion of the function (p(b, mo) becomes more and more difficult since the depen- 
dence of the obsevables on mo is weakening. For example, the time dependent 
magnetization M(t,mo) can be written as M(t, mo) = M^(t,mo) + M(t, 1). The 
part depending on m is only Md(t,m ). For relative big m , e.g. m > 0.50, 
Md(t, mo) is much smaller than M(t, 1). Therefore a small relative error of M(t, 1) 
may induce bigger errors for the function ip(b,mo) when it is determined from 
the scaling plot of M(t,m ). In Fig. |2b], such a situation is shown for a scaling 
plot for L = 144 with m = 0.50. The quality of the scaling collapse is not so 
good. 

The errors of M(t, 1) have many origins. There are apparently statistical 
errors. In the rescaling of M(t,mo), the error of the dynamic critical exponent 
z also enters. More important is the effect of the microscopic time scale t m i C . 
As it is discussed in the numerical measurements of the critical exponent 9 in 
Ref |12|. [16], t m i C becomes bigger as the initial magnetization mo increases. For 
M(t, m = 1), it is observed that t mic ~ 50 JIB] . Finally, the systematic errors of 
the random numbers might also play some role. 

As pointed out in the previous letter, to reduce the extra errors from M(t, 1) 
one can directly estimate the funcion ip(b, mo) from the scaling collapse of M^(t, mo). 
In Fig. ^ a scaling plot of Md(t, mo) for lattice size L = 144 with mo = 0.50 and 
lattice size L = 72 is plotted. Now, the fitting of the data is much better than that 
for M(t, mo). The crosses lie nicely on the fitted solid line in the whole time inter- 
val. As it will carefully be discussed below, the microscopic time scale t m i C shows 
to be quite small here. The resulting value of y?(6, mo) is <p(2, 0.50) = 0.6586(31). 

In any case, errors induced by the error of the exponent z still remain for the 
scaling plot of M^(t,mo). In the above estimation of the values for the function 
(p(b,mo), these errors have not been taken into account. For this reason, we 
have performed the scaling fitting with two different values of z, z = 2.176 and 
z = 2.216. The results for ip(2, 0.50) are (p(2, 0.50) = 0.6571(29) and (p(2, 0.50) = 
0.6602(30) respectively. Compared with the value obtained with z = 2.196, the 
difference is not big. Actually all the values are covered by the statistical errors. 
In contrast to this, the result for the exponent z given in Ref \T2\ is z = 2.196(8) 
and the statistical error is already negligible. 
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In order to reveal the fine struture of the short-time behaviour of the critical 
dynamics, we can proceed in a different way to measure the critical characteristic 
function (p(b, m ). Especially we would like to pay attention to the possible effect 
of the microscopic time scale t mic discussed in Sec. 1, and see also how stable the 
scaling form in Eq.(||) is with respect to the evolution time. Theoretically, the 
scaling form in Eq.([3]) is valid for every time t after the microscopic time scale 
t mic . In Fig. the function ip(b,m ) is determined by searching for the best fit 
of the data in a time interval [10,400]. We may call this 'global' fit. Actually, 
the function ip(b,m ) can be estimated from the data in any 'local' time t. For 
example, we can perform the fitting at a time t in the time interval [t, t + 15] and 
measure the function ip(b,mo). We may call this 'local' fit. 

In Fig. |4a], fit] and |4c], the resulting function (p(2, mo) from the magnetization 
difference Md(t, mo) for the lattice size L = 144 and L = 72 as a function of the 
time t is plotted for the initial magnetization mo = 0.14, 0.30 and 0.60. The errors 
are estimated by dividing the samples with different initial configurations into four 
groups. Taking into account the errors, from the figures we see that y)(2,mo) as 
a function of time t fluctuates around a constant value after a microscopic time 
scale t m i C < 10. This is a strong indication that the scaling form in Eq.(^) is 
valid. Such a small t m t r is consistent with what we observe in the numerical 



measurement of the critical exponent 6 [12|. The final value of cp(2, mo) can 
be estimated by taking the average over a certain time interval. For example, 
for mo = 0.14 we obtain 93(2, 0.14) = 0.1801(6) by averaging over time interval 
[10,400]. This coincides well with that obtained from the 'global' fit. From Fig. 
|4"a| , |4b] and [3c, we can also see the tendency that the fluctuations become bigger 
as mo increases. For m$ = 0.14 the fluctuations are around one percent but for 
m = 0.60 those are already around five percent. 

Here we should stress that the data points of the function ip(2, m ) at dif- 
ferent time t presented in Fig. HaL HE and He] are not completely statistically 
independent. The data points within a time distance of 15 are obtained from 
time intervals which partly overlap each other. Even for the data points with a 
time distance bigger than 15, they are in some way correlated due to the large 
time correlation length. Therefore, using the data in Fig.|4a|, [4b| and f|c| to esti- 
mate the statistical errors is complicated. For these reasons, our statistical errors 
are always estimated by dividing the samples with different initial configurations 
into four groups. However, we should also keep in mind that the data points 
presented in Fig. |4a|, [45] and |4c] are also not completely statistically dependent. 
How these data points correlate in the time direction just represents the charac- 
teristic property of the time evolution of the magnetization. This kind of pictures 
as Fig.|4a|, |4b| and |4c] have shown to be useful in handling the possible effects of 
the microscopic time scale t m i C and of finite lattice sizes |i~9| , Q| . It helps us to 
decide in which time interval we should perform the 'global' fit to obtain the final 
results for the function cp(b, m ). 

The function ip(b,m ) can also be obtain from other observables. In Fig.^j, a 
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scaling plot is shown for the Binder-type cumulant Ud{t,mo) for the lattice size 
L = 144 with m = 0.14 and the lattice size L = 72. The result if (2,0.14) = 
0.1795(8) is as good as that estimated from Md(t, mo). To determine the function 
(p(b, m ) from the cumulant Ud(t, m ), one does not need the static exponent (3/v 
as an input. This is an advantage. The fact that if(b,mo) obtained from both 
Md(t,m ) and Ud{t,m ) coincide well provides strong support for the generalized 
scaling form in Eq.(]3|). However, for bigger initial magnetization, e.g. m > 0.50, 
Ud(t,m ) fluctuates much more than M d (t,m ) does. 

In Table [I], the functions if(2,m ) and if (4, mo) obtained from a global fit 
of the magnetization difference Md(t,mo) are listed. To estimate if (2, mo), the 
fitting is carried out in a time interval [10, 200] for smaller mo, while in a time 
interval [10, 100] for mo > 0.50. For if (4, mo) the fitting is performed in a time 
interval [10, 100] for smaller mo, while in a time interval [10, 50] for mo > 0.50. 
The errors are estimated by dividing the samples with different initial configura- 
tions into four groups. Compared with the data in the previous letter the results 
for bigger m are improved. 

In order to get a more direct understanding of the full critical characteristic 
function if(b,m ), in the previous letter we have defined an effective dimension 
x(b,m ) of the initial magnetization m by if(b, mo) = b x ( fe ' m °) m . Taking the 
results from Table [l], the corresponding effective dimensions x(b,mo) are plotted 
in Fig. ^ Circles are of x(2, mo) and crosses are of x(4, mo). They clearly show 
that the function if(b,m ) is not a constant with respect to m . Later we will 
come back to this point at the end of the next subsection. 

Before closing this subsection, we would like to discuss on the technique of 
the sharp preparation of the initial magnetization in the numerical simulations. 
Taking the Potts model as an example, we randomly generate the initial spin 
configurations such that each spin takes value +1 with probability (1 + 2mo)/3, 
while and —1 with probability (1 — mo)/3. If the lattice size is infinite, a 
sharp value mo for the initial magnetization is automatically achieved. However, 
practically our lattice size is finite. Then the initial magnetization obtained in this 
way fluctuates around the value mo, depending on the lattice size L. This is an 
extra kind of finite size effects. To reduce this effect, we usually adopt the sharp 
preparation technique: when the initial magnetization does not take the value mo, 
we randomly select a spin on the lattice and flip it if the updated magnetization 
comes nearer to mo; we continue this procedure until the initial magnetization 
reaches the value m . In the previous numerical simulations, this technique has 
been widely used and shown its efficiency [TO], [IB], [T8], [TS], |27] . For sufficient big 



lattice sizes, whether one applies the sharp preparation technique does not make 
a big difference, but the data with sharply prepared mo show in general less extra 
finite size effects. How important the sharp preparation technique is depends also 
on the initial magnetization m and what kind of observables one measures. The 
smaller (or nearer to one) the initial magnetization m is, the more important 
sharp preparation technique becomes. In two very recent preprints [14", E0I , more 
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detailed discussions on this problem have been presented. A typical example is 
the measurement of the global persistence exponent |23], |20|. Without the sharp 
preparation technique, it is very difficult to observe the power law decay of the 
persistence probability of the (global) magnetization and the extra finite size 
effect remains even when the lattice size L = 256 for the two dimensional Ising 



model 120 



For the measurement of the critical characteristic function, we have also re- 
peated some calculations without the sharp preparation of the initial magnetiza- 
tion. In Fig. 0, a scaling plot for the magnetization difference M d {t) is displayed 
for the lattice size L = 72 and L = 144. The initial magnetization for L = 144 is 
mo = 0.14. Stars represent the time evolution of the magnetization difference for 
L = 144, while the lower and upper solid line are those for L = 72 with mo = 0.175 
and mo = 0.185 respectively. For comparison, the correponding magnetizations 
for the lattice size L = 144 and L = 72 with the sharp preparation of the initial 
magnetization have also been displayed by the dotted line and the dashed lines 
in the figure. We see the difference between with and without the sharp prepara- 
tion technique is already not big. The crosses are the same data as the stars but 
rescaled in time. The resulting value of tp(b, mo) is tp(2, 0.14) = 0.1811(9). Within 
the statistical errors, this is consistent with (^(2,0.14) = 0.1800(4) obtained with 
the sharp preparation technique. Similarly, without the sharp preparation tech- 
nique we obtain <p(2,QA0) = 0.5283(17). Looking at Table |, all these results 
show that for the lattice as big as L = 72 and L = 144 the difference between 
with and without the sharp preparation technique is within the statistical errors. 

3.2 The Ising model 

The Hamiltonian for the Ising model is 

H = JY,Si Sj , Si = ±l, (10) 

<ij> 

with < ij > representing nearest neighbours. In equilibrium the Ising model 
is exactly solvable. The critical point locates at J c = log(l + a/2)/2 and the 
exponent (3/v = 1/8. The kth moment is defined as 

M(t) = (l^(T,s l (t)] k ). (11) 

i 

In this section, we follow the same procedure as demonstrated in the last 
subsection to study the short-time behaviour of the critical dynamics for the 
Ising model and to determine the critical characteristic function for the initial 
magnetization. It has been pointed out in the previous letter [p5jl , for the Potts 
model in the determination of if (2, mo) the finite size effect for the lattice pair 
(Li, L 2 ) = (72, 144) is already very small. This situation is similar for the Ising 
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model. Therefore, in this section we restrict our discussions to the lattice pair 
(Li,L 2 ) = (64,128). For the determination of <£>(4, m ), we use the lattice pair 
(Li, L 2 ) = (32, 128). For lattice sizes L = 36, 72 and 144, we have performed sim- 
ulations with 80 000, 32 000 and 8 000 samples of different initial configurations 
respectively. The initial magnetization for the lattice L = 128 takes the value 
m = 0.10, 0.20, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80 and 1.00. For the lattice 
L = 128 with m > 0.60, the statistics is 16 000. The maxmum evolution time 
is t = 2500 for the lattice size L = 128, t = 250 for the lattice size L = 72 and 
t = 100 for the lattice size L = 36. 

In case of the Potts model, we believe that the dynamic exponent z = 2.196 
used in the last subsection is very close to the true value [[L(J [12], |j~6| . However, 



for the Ising model the value of the dynamic exponent is still a matter of debate. 



Values ranging from z = 2.13 to z = 2.17 can be found in the literature [28, 29 



50] , [18], [19|, |3T], |32| . Fortunately, in the determination of the critical characteristic 
function <p(b, mo) the deviation about one percent in the value of z does not induce 
a big change in the results for (p(b,m ), at least for not too big m . In Fig. || 
a scaling plot for the magnetization of the lattice size L = 128 with itlq = 0.40 
and lattice size L = 64 is displayed. Crosses and circles are the same data for the 
magnetization of L = 128, but rescaled with two different values of the dynamic 
exponent z, z — 2.15 and z = 2.17 respectively. We can see the difference is very 
small. The resulting values of <£>(2,mo) are <£>(2,0.40) = 0.5588(25) for z = 2.15 
and (p(2,0A0) = 0.5567(25) for z = 2.17. On the other hand, in our numerical 
measurements from the short-time dynamics the value of z tends to z = 2.15 
[ |19| , [i~2|j . Therefore, in this section we always use the value z = 2.15 in the 
analyses of the data unless we specify other values. Later we will give some more 
details on this point. 

Here we should also mention that as in the case of the Potts model, for the 
lattice pair (Li,L 2 ) = (64,128) the difference between data with and without 
the sharp preparation technique, i.e. the extra finite size effect from the initial 
configurations is already not big. For example, without the sharp preparation 
of the initial magnetization we get <^(2,0.40) = 0.5566(18) from the scaling col- 
laps of the magnetization. Wihtin the statistical errors, this is consistent with 
<f(2, 0.40) = 0.5588(25) obtained with the sharp preparation technique. 

In Fig. H a scaling plot for the magnetization of the lattice size L = 128 with 
mo = 0.30 and lattice size L = 32 is displayed. The scaling collapse looks also 
very nice even though the fluctuation is slightly bigger than that in the scaling 
plot for the lattice pair (L 1; L 2 ) = (64,128). The resulting value of <£>(2,m ) is 
<f (4, 0.30) = 0.585(2). Practically, the measurement of y?(4, m ) is more difficult 
than that of ip(2, mo). We have to perform simulations for bigger lattice sizes and 
longer evolution time. Much computer time is needed. Therefore, to measure 
ip(b, m ) for even bigger b, it is not wise to follow the way discussed in this paper. 

To see the possible effect of the microscopic time scale t mic , we can analyse 
the data also by the local fit. In Fig. [1C], ip(2, m ) with m = 0.40 obtained from 
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the magnetization with lattice size L = 128 and L = 64 is plotted versus time 
t. We can see clearly that t mic is also very small, t mic < 10. Therefore, in this 
section we always start our measurements of the critical characteristic function 
<f (2, mo) from the time t = 10. 

In case of the Potts model, when mo is bigger the scaling collapse for the 
magnetization M(t, mo) is not as good as that for the magnetization difference 
Ma(t,mo). However, it is interesting that for the Ising model both scaling plots 
for M(t, m Q ) and Md(t, mo) give similar results, at least when one takes z = 2.15. 
This shows that in the simulations we do not get so much fluctuations from 
M(t, 1). In other words, the data of M(t, 1) for the lattice pair (Li, L2) collapse 
very well. In Table [| the function ip(2, mo) and ip(4, mo) obtained from a global 
fit of the magnetization M(t,m ) are listed. For z = 2.17, the case is slightly 
different. We will discuss this later. 

In Fig. |6], the final results for the effective dimension x(2, mo) and x(4, mo) 
are given. The data are taken from Table |2|. Stars with solid line are of x(2, mo) 
and diamonds with dashed line are of x(4, mo). From Fig. ^ we observe that 
the effective dimension x(b, mo) for the Ising model and the Potts model is quite 
different. We would expect that x(b,m ) always monotonously decreases as m 
goes from m = to m = 1. This is indeed the case for the Ising model. 
However, for the Potts model x(b, mo) first increases and after some time reaches 
its maximum and then decreases. This is consistent with what one observes in 



the numerical measurement of the exponent 9 [12 . When mo increases from 



mo = 0.0, for the Ising model the effective 9 decreases while for the Potts model 
it increases. This phenomanon remains to be understood. It is challenging to 
derive analytically the generalized dynamic scaling form in Eq.(|3[) and the critical 
characteristic function ip(b, mo). 

Now, naturally one may ask, is the critical characteristic function ip(b, mo) 
really non-trivial? In other words, it may happen that if we take another field 
rather than the initial magnetization as the scaling variable one exponent is suffi- 
cient to describe the scaling behaviour of the initial magnetization. Theoretically, 
this is possible. For example, we can take g = F _1 (m-o) as a scaling variable as 
it is discussed in Sec. 2. However, the problem is not so simple. First, in general 
we don't know the function F~ 1 (mo)- Secondly, the physical meaning of g is 
not clear. Actually, F _1 (mo) is just the function we intend to measure. Very 
recently, we have received a preprint [E6J by U. Ritschel and P. Czerner (Univ. 
Essen, FRG) where a comment is given to the topic addressed in this paper and 
the previous letter ||25|| . It is said that if one uses the initial external magnetic 
field as the scaling variable rather than the initial magnetization, the exponent 
Xo is sufficient to describe the dynamic scaling. They got an indication for that 
from their data of x(b,mo) for the Ising model. However, their conjecture ap- 
parently fails in the Potts model. For the Ising model, their data overlap partly 
and are consistent with ours within the errors. But our lattice sizes are bigger 
and our errors are smaller. Our data within errors are not consistent with their 
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conjecture. Actually, some points of their data are also outside of their expected 
values. 

At high temperature, the external magnetic field h and the initial magnetiza- 



tion mo have a simple relation, for the Ising model |26] 

777 = tanh(/i), (12) 
and for the Potts model ^ 

mo = 3 2 i ^-l ) - (13) 
Supposing h is used as a scaling variable in the scaling form in Eq.(|J), the effective 
anomalous dimension y(b, mo) of h can also be extracted from the numerical data. 
In the limit 777 = 0.0, y(b,0) = x(b,0) = Xq. The exponent xq, which has been 
measured in Ref. |I| [12| , is x = 0.536 and x = 0.298 for both the Ising and 
the Potts model respectively. The results are shown in Fig. [TT| For the Potts 
model, 7/(6, mo) apparently is not a constant. For the Ising model, it is not too 
far from a constant but there is significant difference. Especially, as is pointed 
out before, the results obtained from both Md(t, mo) and M(t,m ) agree very 
well. For mo > 0.80, the fluctuations for the Potts model become very big. We 
can not see whether 7/(2, mo) continues to increase or starts to decrease. For the 
Ising model, the data are not available. 

We have discussed that the dynamic exponent z is not known so rigorously. 
Therefore we may also present the results with other values of z. For example, 
if one takes z = 2.17 for the Ising model, the resulting 7/(2, m ) is displayed in 
Fig. [l!| Stars are of 7/(2, mo) measured from M(t, mo) and squares are that from 
Md(t,mo). The results obtained from Md(t,m ) coincide very well with those 
measured with z = 2.15 (the solid line). However, for relative bigger m the 
data from M(t,mo) show some deviation. We believe the result of M^(t,mo) are 
nearer to the true values of 7/(2, mo). In the scaling plot for M(t,mo) one may 
have got some errors from M(t, 1). This might also indicate that z = 2.15 is a 
better value for the dynamic exponent z of the Ising model. Anyhow, even the 
data from M(t,m ) are not compatable with a constant. 



4 Conclusions and remarks 

We have numerically simulated the universal short-time behaviour of the dynamic 
relaxation process for the two dimensional Ising and the Potts model at criticality 
starting from an initial state with very high temperature and arbitrary magneti- 
zation. The critical characteristic function for the initial magnetization, which is 
introduced to describe the generalized dynamic scaling, is determined with finite 
size scaling. It is interesting that the critical characteristic functions for the Ising 
and the Potts model look quite different. When the initial magnetization varies 
from m = to m = 1, the effective dimension of the initial magnetization 
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for the Ising model monotonously decreases, while that for the Potts model first 
rises and then decreases. It is challenging to derive analytically the generalized 
scaling form as well as the critical characteristic function ip(b,m ). However, the 
standard perturbative renormalization group methods for the continuum models, 
as e.g. the 4 -theory, do not apply here. 

The initial magnetization can be considered as a coupling of the systems. The 
generalized scaling form describes the crossover scaling behaviour of the dynamic 
systems between the fixed point m = to m = 1. This kind of phenomena 
can be found in many statistical systems, for example, at the ordinary phase 
transition in the surface critical phenomena. When a non-zero magnetization is 
applied to the surface, a similar scenario as in the short-time dynamics appears. 

Acknowledgement: The author would like to thank L. Schiilke for helpful 
discussions. 
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too 




X{%m ) 


n 1 /i 
U.14 


U.loUU(U4J 


U.zoo(l ) 


0.22 


0.2866(10) 


0.378(2) 


0.30 


0.3952(11) 


0.521(3) 


0.40 


0.5299(07) 


0.691(4) 


0.50 


0.6586(31) 


0.815(8) 


0.60 


0.7678(23) 


0.907(6) 


0.70 


0.8657(29) 


0.954(3) 


0.80 


0.9272(91) 


0.987(7) 



Table 1: The characteristic function tp{b, too) for the Potts model measured from 
the magnetization difference Md(t). 



too 


y?(2,TO ) 


</?(4,m ) 


0.10 


0.1439(04) 


0.206(01) 


0.20 


0.2857(14) 


0.405(03) 


0.30 


0.4254(10) 


0.585(02) 


0.40 


0.5588(25) 


0.741(07) 


0.50 


0.6800(31) 


0.847(13) 


0.60 


0.7890(38) 


0.931(13) 


0.70 


0.8783(78) 


0.964(12) 



Table 2: The characteristic function (p(b, m ) for the Ising model measured from 
the magnetization. 
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Figure 1: Magnetization profiles for the Potts model with L = 72 and tuq = 
0.14, 0.22, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80, 0.90 and 1.00. 
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Figure 2a: The scaling plot for the magnetization with (Li,L 2 ) = (72,144) 
for the Potts model. For the lattice size L 2 = 144, the initial magnetization 
m = 0.14. 
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Figure 2b: The scaling plot for the magnetization with (Li,L 2 ) = (72,144) 
for the Potts model. For the lattice size L 2 = 144, the initial magnetization 
m = 0.50. 
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Figure 3: The scaling plot for the magnetization difference M d (t) with (Li, L 2 ) = 
(72, 144) for the Potts model. For the lattice size L 2 = 144, the initial magneti- 
zation m = 0.50. 
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Figure 4a: The characteristic function </?(2,m ) estimated with (Li,L 2 ) = 
(72, 144) from a 'local' fit of M d {t) for the Potts model. For the lattice size 
L 2 = 144, the initial magnetization m = 0.14. 
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Figure 4b: The characteristic function </?(2,m ) estimated with (L 1 ,L 2 ) = 
(72, 144) from a 'local' fit of M d {t) for the Potts model. For the lattice size 
L 2 = 144, the initial magnetization m = 0.30. 
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Figure 4c: The characteristic function (p(2,m ) estimated with (L 1? L 2 ) = 
(72, 144) from a 'local' fit of M d {t) for the Potts model. For the lattice size 
L 2 = 144, the initial magnetization m = 0.60. 
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Figure 5: The scaling plot for the Binder-type cumulant U d (t,m ) with 
(L U L 2 ) = (72,144) for the Potts model. For the lattice size L 2 = 144, the 
initial magnetization m = 0.14. 
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Figure 6: The effective dimension x(b,mo) for both the Ising and the Potts 
model. For the Potts model the data are obtained from the magnetization dif- 
ference M d (t,m ) with (Li,L 2 ) = (72, 144) and (36, 144). Circles with solid line 
are of x(2,m ) and crosses with dashed line are of x(A, mo). For the Ising model 
results are obtained from the magnetization M(t,mo) with (L\,Li2) = (64, 128) 
and (32,128). Stars with solid line are of x(2, m ) and diamonds with dashed 
line are of x(4, m ). 
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Figure 7: The scaling plot for the magnetization difference M d (t) with (L 1 , L 2 ) = 
(72, f44) for the Potts model. For the lattice size L 2 = 144, the initial mag- 
netization m = 0.14. The sharp preparation technique is not adopted. For 
comparison, the corresponding magnetization difference profiles with the sharp 
preparation technique are also displayed for the lattice sizes L = 144 and L = 72 
with the dotted line and the dashed lines respectively. 
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Figure 8: The scaling plot for the magnetization with (Li,!^) = (64,128) for 
the Ising model. Crosses correspond to z — 2.15 and circles are from z = 2.17. 
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Figure 10: The characteristic function </?(2,m ) estimated with (Li,L 2 ) = 
(64, 128) from a 'local' fit of M(t) for the Ising model. For the lattice size 
L 2 = 128, the initial magnetization m = 0.40. 
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Figure 11: The effective dimension y(2,m ) for both the Ising and the Potts 
model. The long dashed and dashed line are the exponent Xq for the Ising and the 
Potts model respectively, measured in the limit m = 0.0. Crosses are of y(2, m ) 
for the Ising model calculated from Table 2 and circles are the same but obtained 
from the scaling plot of M d (t,m ). Triangles are of y(2,m ) for the Potts model 
obtained from Table 1. 
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Figure 12: The effective dimension 7/(2, rao) for the Ising model with z = 2.17. 
The long dashed line is the exponent x for the Ising model. Stars are of 7/(2, m ) 
calculated from M(t, m ) and the squares are the same but obtained from the 
scaling plot of M d (t,m ). The solid line is that estimated from M d (t,m ) but 
with z = 2.15. 
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